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Abstract 

We consider compact metric graphs with an arbitrary self adjoint realisation of the 
differential Laplacian. After discussing spectral properties of Laplacians, we prove 
several versions of trace formulae, relating Laplace spectra to sums over periodic 
orbits on the graph. This includes trace formulae with, respectively, absolutely and 
conditionally convergent periodic orbit sums; the convergence depending on proper- 
ties of the test functions used. We also prove a trace formula for the heat kernel and 
provide small-t asymptotics for the trace of the heat kernel. 
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1 Introduction 



Some ten years ago Kottos and Smilansky |KS97l IKS99bj introduced quantum graphs as 
convenient models in the field of quantum chaos [HaaOll IStoQQ] . where a major goal is 
to understand the connection between dynamical properties of a quantum system and 
its associated classical counterpart |BGS84j . Previously introduced models that possess 
classical counterparts with chaotic dynamics include quantum billiards, motions on Rie- 
mannian manifolds with negative curvatures and quantum maps. These models have been 
studied with considerable success, however, they have often turned out to bear unwanted 
complications. Quantum graphs are constructed along the lines of many of these mod- 
els in that they are mainly concerned with spectral properties of Laplacians. In a sense 
they are maximally reduced versions of such models in that the underlying configuration 
space is one dimensional. The non trivial topology of the graph, however, introduces suffi- 
cient complexity such that the quantum system behaves like typical quantum systems with 
chaotic classical counterparts, see |KS97t IKS99b] . On the other hand, many details of the 
classical dynamics are considerably simpler and quantum spectra are generally known to 
some more detail, so that quantum graph models proved to be very useful in the field of 
quantum chaos |GS06] . 

Trace formulae provide a direct connection between classical and quantum dynamics 
in that they relate quantum spectra to classical periodic orbits. In general, this connec- 
tion arises in the form of an asymptotic relation, valid for large wave numbers. There 
are only few exceptional cases in which trace formulae are identities. Among these are, 
most notably, Laplacians on flat tori and on Riemannian manifolds with constant negative 
sectional curvatures. In both cases the configuration manifolds are Riemannian symmetric 
spaces, allowing for the application of powerful methods of harmonic analysis in proving the 
relevant trace formulae, i.e., the Poisson summation formula and the Selberg trace formula 
|Sel56] . respectively. Generically, however, the tools of harmonic analysis are not avail- 
able, and trace formulae have to be proven using semiclassical or microlocal techniques, 
which naturally involve asymptotic methods. Semiclassical trace formulae were introduced 
by Gutzwiller jGutTlj for the spectral density of quantum Hamiltonians. Subsequently, 
Balian and Bloch [BB72j set up analogous trace formulae for Laplacians on domains in 
M" in a short wavelength approximation. The first mathematical proofs, for Laplacians on 
Riemannian manifolds, are due to Colin de Verdiere |CdV73j as well as Duistermaat and 
Guillemin |DG75j . Later, proofs for the semiclassical case followed |Mei92t rPU95j . 

One of the virtues of quantum graph models that led Kottos and Smilansky to intro- 
duce them to the field of quantum chaos is that their trace formulae are identities, very 
much in analogy to the Selberg trace formula. The first quantum graph trace formula, 
however, is due to Roth |Rot83j . who expressed the trace of the heat kernel for a Laplacian 
with Kirchhoff boundary conditions in the vertices of the graph as a sum over periodic 
orbits. Kottos and Smilansky |KS99bj then introduced a trace formula for the spectral 
density of the Laplacian, and later Kostrykin, Potthoff and Schrader [KPS07] extended 
Roth's trace formula to more general boundary conditions. In these cases the boundary 
conditions characterising the domain of the Laplacian were of a non-Robin type in that 
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they do not mix boundary values of functions and their derivatives. This leads to periodic 
orbit sums in the trace formulae that closely resemble those occurring in the Selberg trace 
formula. In this paper our principal goal now is to consider general self adjoint realisations 
of Laplacians on compact metric graphs and to prove associated trace formulae with fairly 
general test functions. In particular, we do not require their Fourier transforms to be 
compactly supported. Allowing for Robin-type boundary conditions leads to trace formu- 
lae that are still identities, yet the amplitudes multiplying the oscillating factors in each 
periodic orbit contribution depend on the wave number in a non trivial way. In principle, 
these amplitude functions are known and possess asymptotic expansions for large wave 
numbers. Therefore, in a sense these trace formulae are intermediate between Selberg 
and Gutzwiller/Duistermaat-Guillemin trace formulae, where in the latter case only the 
asymptotic expansions of the amplitude functions are known, and the test functions must 
have compactly supported Fourier transforms. 

This paper is organised as follows: In section [2] we briefly review the construction of 
quantum graphs, including parametrisations of self adjoint realisations of the Laplacian 
developed by Kostrykin and Schrader |KS99a] as well as Kuchment jKuc04j . Following 
this we investigate various properties of quantum graph edge S-matrices as introduced 
in |KS99bj . focussing on their analytic properties. In section H] we then discuss general 
properties of Laplace spectra on compact graphs. The main part of this paper can be found 
in section [5] where we prove several versions of quantum graph trace formulae: Theorem l5.3l 
contains a trace formula with a double sum over periodic orbits, of which only the sum 
over topological orbit lengths converges absolutely, and which allows for a large class of 
test functions. In Theorem 15.41 however, we restrict the class of test functions with the 
effect that the entire sum over periodic orbits converges absolutely. A suitable choice of a 
test function then allows to establish the trace of the heat kernel for arbitrary self adjoint 
realisations of the Laplacian, see Theorem 15.51 We finally summarise and discuss our 
results in section [Ul Some of the results presented here have been announced in jBEOSj . 

2 Preliminaries 

We begin with reviewing the relevant concepts underlying the construction of quantum 
graphs. 

2.1 Metric graphs 

In the sequel we shall consider finite, metric graphs F = (V, I). Here V is a finite set of 
vertices {f i, . . . , wv} and £ is a. finite set of edges {ei, . . . , e^;}. When an edge e connects 
the vertices v and w, these are called edge ends. Two edges are adjacent, if they share an 
edge end; loops and pairs of multiply connected vertices shall be allowed. The degree of 
a vertex v specifies the number of edges with v as one of their edge ends. A metric structure 
can be introduced by assigning intervals [0, /j] to edges Cj, along with coordinates Xi G [0, U]. 
The i?-tuple I = {li, . . . ,Ie) then collects all edge lengths. Through this procedure the edge 
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ends are mapped to the end points of the intervals in a specified manner. The edge end 
corresponding to Xj = is then called initial point and, correspondingly, the other end 
point is the terminal point of the edge Cj. Hence not only the connectedness of the graph 
is specified, but also an orientation of the edges is introduced. We emphasise that the 
specific choice of the orientation thus made does not impact the results of this paper, in 
fact the choice of the initial and terminal points of the edges are arbitrary. The intervals 
are only used to construct the Laplace operator as a differential operator. 

It is useful to arrange the 2E edge ends in a particular way: we list the initial points 
in the order as they occur in the list of edges, followed by the terminal points in the same 
order. For the trace formula we, moreover, require the following notions. 

Definition 2.1. A closed path in F is a finite sequence of edges (ej)^^j^, such that the edges 
Cj and Cj+i for G {1, . . . ,n — 1} and the edges and ei are adjacent. A periodic orbit is 
an equivalence class of closed paths modulo cyclic permutations of the edges. 

The number n of edges in a periodic orbit p is its topological length, whereas the sum 
Ip = hi + ■ ■ ■ + hn of the metric lengths of its edges is the metric length, or simply length, 
of p. A periodic orbit is primitive, if it is not a multiple repetition of another periodic 
orbit. Furthermore, the set of periodic orbits of the graph is called V, and P„ is its subset 
of orbits with topological length n. 



2.2 Quantum graphs 

Quantum mechanics on a metric graph can be studied in terms of the Schrodinger equation 

i^^, = -A^, . (2.1) 

Here ipt is a vector in an appropriate Hilbert space and A is a suitable (differential) Laplace 
operator, acting on functions defined on the edges of the graph. Since this operator is of 
second order, it is independent of the choice of the coordinates Xj or Xj := Ij — xj on the 
edges Cj, and therefore also independent of the specification of edge ends as initial and 
terminal. 

In order to be in a position to view the negative Laplacian in (12.11) as a quantum 
Hamiltonian, it has to be realised as a self adjoint operator on a suitable Hilbert space. 
We therefore introduce the space of square integrable functions on F as the quantum graph 
Hilbert space. In this context a function on the graph is a collection F = (/i, . . . , /e) of 
functions fj : (0, Ij) — C on the edges. Therefore, one considers the following function 
spaces, 

E E 

C-(F) = 0C°°(O,/,) and L\r) = ^ L\0,l,) . 

The latter is a closed orthogonal sum of Hilbert spaces with respect to the scalar product 

^ rh 

{F,G) -.= ^2 / fji^j)9jixj) dxj . 
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As a differential expression the (negative) Laplacian is simply given by 

-AF:=(-/f,...,-/^), 

where dashes denote derivatives. This expression may now serve as a way to introduce a 
closed, sjTiimetric operator {—A,Vq) with domain 

E 

i=i 

Here each term in the orthogonal sum consists of an L^-Sobolev space of functions which, 
together with their derivatives, vanish at the edge ends. The deficiency indices of this 
operator are {2E,2E), and thus it possesses self adjoint extensions that can be classi- 
fied according to von Neumann's theory (see, e.g., |RS75j ). An alternative approach has 
been developed in detail by Kostrykin and Schrader |KS99a] . which provides a convenient 
parametrisation that is particularly useful for later purposes. In this context one introduces 
the boundary values 

F,, = (/i(0), . . . , fE{0), flih), . . . , fEilE)f , 

FL = (/((o), . . . , /^(o), -f[{h), . . . , -fE{lE)f , 

of functions and their derivatives, whereby the signs ensure that inward derivatives are 
considered at all edge ends. Notice that the order of the terms follows the convention 
of arranging edge ends introduced previously. Boundary conditions on the functions in 
the domain of a given self adjoint operator are specified through a linear relation between 
boundary values; they are of the form 

AFb, + BFi = , (2.3) 

see |KS99aj . Here A,B & M{2E, C) are two matrices such that 

• the matrix [A, B), consisting of the columns of A and B, has maximal rank 2E, 

• AB* is self adjoint. 

These conditions then imply the self adjointness of the operator, and every self adjoint 
extension can be achieved in this manner. Occasionally, we shall denote a particular such 
self adjoint realisation of the Laplacian as A{A,B,l). This parametrisation, however, is 
obviously not unique because a multiplication of (12.31) with C e GL(2ii^,C) from the left 
does not change the boundary conditions. On the other hand, if A[A, B, I) = A{A', B', I), 
there exists C G GL(2E,C) with A' = CA and B' = CB, see |KS99a] . Thus, for any 
C e GL(2ii^,C) both A,B and A' = CA,B' = CB provide an equivalent characterisation 
of the same operator. 

The linear relations (12.31) can in principle relate boundary values at any set of edge ends. 
We wish, however, the operator to respect the connectedness of the graph and therefore we 
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restrict ourselves to local boundary conditions. These are characterised by the condition 
that (12.31) only relates edge ends that form a single vertex. To this end we now group the 
edge ends in (12. 2p according to the vertices they belong to. Local boundary conditions 
then lead to a block structure of the matrices A and B, 



A = A, and 5 = 05,, 



(2.4) 



such that each block, represented by A^ and B^, exactly relates the boundary values of 
functions and their derivatives at the vertex v. In this context self adjointness of the 
Laplacian is achieved, if for all f G V the rank of {A^, B^) is and A^B* is self adjoint. 

To mention a few examples, a vertex with Dirichlet boundary conditions can be char- 
acterised hj A^ = By = 0, whereas for Neumann boundary conditions one would 
choose Ay = 0, By = 1^^. Moreover, the generalised Kirchhoff boundary conditions used 
by Kottos and Smilansky |KS99b] can be achieved by choosing 



A. 



/I 



v 



-1 



and By 



( 



\ 



V 



(2.5) 



where only the non vanishing matrix entries are indicated. Here [ly must be real; when 
/i„ = the usual Kirchhoff conditions are realised. 

The approach we have taken here when singling out local boundary conditions is to 
start from a given graph and then to realise the Laplacian as a self adjoint operator on that 
graph. An alternative view would be to consider the above construction as a realisation of 
the Laplacian as a self adjoint operator on a collection of E intervals. One can then find a 
graph, i.e., a way to connect the interval ends in vertices, such that a given pair of matrices 
A and B is local with respect to this graph; there even is a unique graph maximising the 
number of vertices [KSOGaj . In the following we shall, however, stick to the previous view. 

The non-uniqueness in the choice of the matrices A and B can be overcome by parametris- 
ing the self adjoint realisations of the Laplacian in terms of projectors onto subspaces of 
the 2£'-dimensional spaces of boundary values. To this end Kuchment |Kuc04j introduced 
the projector P onto the kernel of B as well as the projector Q = \—P onto the orthogonal 
complement (ker i?)-*- = ran B* in C^^ and proved that A maps ran B* into ran B. He then 
defined the (self adjoint) endomorphism 



'■— (-B|ranB*) AQ 

of rani?*, and showed that the boundary conditions (12. 3p are equivalent to 



bv 







and LQFby + QF^. 



bv 







(2.6) 



(2.7) 



Moreover, there exists a C G GL{2E, C) such that 

A' = CA = P + L and B' = CB = Q, 
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implying that 



L = A'B' 



A refinement of this construction can be found in |FKW07] . From fl2.7p one concludes that 
in cases where L = 0, the boundary conditions do not mix boundary values of the functions 
themselves with those of their derivatives. We call these non-Robin boundary conditions, 
and all other cases Robin boundary conditions. 

Yet another description of the self adjoint realisations of the Laplacian employs the 
associated quadratic forms. As shown by Kuchment |Kuc04j . given a realisation of the 
Laplacian as a self adjoint operator, the associated quadratic form 

Qa[F] = V/ \f',{x)\^ dx-{F,,,LF!^)c^E (2.9) 

has a domain that consist of all functions F on the graph, with components /e G -f^^(0, Ig), 
whose boundary values fulfil PF^v = 0. 

3 The S-matrix 

In quantum graph models Laplace eigenvalues can be conveniently characterised in terms 
of zeros of finite dimensional determinants, and thus these models are amenable to pow- 
erful analytical as well as numerical methods. In quantum billiards a related method was 
pioneered by Doron and Smilansky |DS92j as the scattering approach to quantisation. In 
general, this method relies on semiclassical approximations. As first demonstrated by Kot- 
tos and Smilansky |KS99bj . however, in quantum graphs the scattering approach allows to 
determine Laplace eigenvalues exactly from a finite dimensional secular equation. 

The scattering approach bears its name from the fact that it is based on scattering 
processes occurring when one opens up a given closed quantum system appropriately. In 
quantum graphs the procedure of opening up consists of replacing each vertex and its 
attached edges by an infinite star graph. This is the single, given vertex v with infinite 
half lines attached that replace the edges of finite lengths. Carrying over the local boundary 
conditions at the vertices from the original closed quantum graph, one thus obtains V open 
quantum systems, which each possess an on-shell scattering matrix a^{k). More precisely, 
this S-matrix is defined in terms of d^ functions F^^^ on the infinite star graph associated 
with V, whose components on the dy infinite edges are 

g-ifex _j_ (j''..{k) e'''^ , j = i 

One then requires each of these functions to fulfil the boundary conditions aX v. In terms 
of the parameterisation of the boundary conditions described in Section 12.21 one then finds 
that 

a^'ik) = -{Ay + ikBy)-\Ay - ikBy) , for A; e M\{0} . (3.1) 
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= -^^^ + Tirr— ■ (3-2) 



The conditions imposed on A^j and By in order to achieve self adjoint boundary conditions 
ensure that ± ikB^ are invertible and that the vertex S-matrix a^{k) is unitary for all 
k G M \ {0}, see |KS99a] . 

The local scattering matrix associated with a vertex with Dirichlet or Neumann bound- 
ary conditions is a"" = — ld„ or = respectively. In contrast, according to fl2.5p 
generalised Kirchhoff boundary conditions lead to a vertex S-matrix with elements of the 
form |KS99bj 

2k 

The local S-matrices of the entire graph can now be grouped together vertex by vertex. 
In that process all edges occur twice, namely associated with the vertex S-matrices of their 
two edge ends. It is hence useful to consider directed edges, and view the matrix elements 
of as describing transitions from a directed edge with terminal point f to a directed 
edge that has v as its initial point. As the list of all directed edges corresponds to listing 
their initial points, the transitions prescribed by can also be performed on the boundary 
values fl2.2p . As a result one obtains the matrix 

S{A, B; k) = ~{A + ikB)-\A - ikB) . (3.3) 

Again, A±ikB are invertible and S{A, B; k) is unitary for all G M \ {0}. Moreover, (13. 3p 
is invariant under the substitution of A, B by CA, CB for all C G GL(2i?, C) and therefore 
is associated with the self adjoint realisations of the Laplacian [KSOGaj . 

We emphasise that, although (13.31) can be defined for any self adjoint realisation, the 
vertex S-matrices (13. ip can only be recovered from (13.31) in the case of local boundary 
conditions, i.e., when (12. 4p is satisfied. Despite the fact that a closed quantum graph does 
not allow for quantum scattering in the usual sense, the quantity (13. 3p is often referred to 
as the edge (or bond) S-matrix |KS99bj of the quantum graph. 

Furthermore, using the parametrisation (12. 7p of boundary conditions and utilising (12. 8p 
as well as the fact that (L + ik)~^ commutes with L — ik, one obtains the representation 

S{A,B;k) = -P -Q{L + ik)-\L -ik)Q , G M\{0} , (3.4) 

for the edge S-matrix, see [KSOGa] . 

From the expressions (13.30 and (13. 4p it appears that the S-matrix generally depends 
on the wave number in a non-trivial way. However, certain boundary conditions lead 
to fc-independent S-matrices. Obvious examples are Dirichlet and Neumann boundary 
conditions, as well as the usual Kirchhoff boundary conditions, i.e., (13. 2p with /i„ = 0. A 
general characterisation of such boundary conditions was provided in |KPS07j in terms of 
the following equivalent conditions: 

• S{A, B; k) is ^-independent. 

• S{A, B; k) is self adjoint for some, and hence for all, /c > 0. 

• S{A^ B; k) = 1 — 2P for some, and hence for all, k > 0. 
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• AB* = 0, i.e., L = 0. 



The last point shows that fc-independent S-matrices arise exactly in the case of non-Robin 
boundary conditions. 

Below we are going to prove some properties of S-matrices that are relevant for the 
trace formula. In this context, for Robin boundary conditions an important role will be 
played by the spectrum cr(L) of the self adjoint matrix L (12. 6p . 

We shall make extensive use of the S-matrix extended to complex wave numbers k and 
therefore need the following. 

Lemma 3.1. Let A and B specify self adjoint boundary conditions for the Laplacian on 
the graph. Then the S-matrix ^3. 3^ has the following properties: 

1. ^(A, B] k) can be continued into the complex k-plane as a meromorphic function, and 
has simple poles at the points of the set icr(L) \ {0}. 

2. S{A, B; k) is unitary for all A; G M. 

3. 5'(A, B] k) is invertible for all k ^ C \ [±icr(L) \ {0}], and its inverse is S{A, B; —k). 

Proof. We henceforth extend the self adjoint endomorphism L of rani?*, see (12.61) . to an 
endomorphism of C"^^ by setting it to zero on (rani?*)-'- = kevB. We then diagonalise 
L utilising an appropriate unitary W, and denote the non-zero eigenvalues (counted with 
their multiplicities) by {Xi, . . . , X^}. This leaves the eigenvalue zero with a multiplicity 
of 2E — d. There are r := dim ran 5* — d (ortho normal) eigenvectors of L in rani?* and 
s := dimker i? = 2E — dim rani?* eigenvectors in ker i?, respectively, corresponding to the 
eigenvalue zero. 

Employing this diagonalisation in the representation (13. 4p of the S-matrix then leads 
to the expression 

/ Ai-ifc \ 
' Ai+ifc » 



S{A,B-k) = W* 



Ad+ifc 



W . (3.5) 



Since the unitary W is independent of k the first statement of the lemma is obvious. 

The unitarity of S{A, B; k) for real k also follows immediately by observing that the 
diagonal entries in (13. 5p are all of unit absolute value. 

The third statement follows in a completely analogous fashion from the representation 
(1331). □ 



Knowing that the S-matrix is analytic in k, one would like to calculate its derivative. 
This in fact is required in the proof of the trace formula below. It is even possible to relate 
the derivative of S{k) to the S-matrix itself. 



9 



Lemma 3.2. Under the same assumptions as in Lemma \3.1\ one obtains for k E C \ 
[±ia(L)\{0}], 

^SiA, B- k) = [S{A, B; k) - SiA, B, k)-'] B; k) . (3.6) 

We remark that for real k the unitarity of the S-matrix can be invoked to obtain from 
(13. 6p that it is independent of fc, iff it is self adjoint. 

Proof. Let us first assume that k G M\ {0} and abbreviate ^(A, B\ k) as S{k). We also de- 
note derivatives w.r.t. by a dash and use the relation ^[X(A;)]^"'^ = —X(k)^^X'(k)X(k)~^, 
which is true for any differentiable function X{k) taking values in GL{2E, C). Recall that 
the conditions imposed on A, B ensure that A ± ikB is invertible for G M \ {0}. Hence 

S'{k) = {A + ikB)-\iB){A + ikB)-\A - ikB) + {A + ikBy\iB) 

= -i{A + ikB)~^B{S{k) - 1) . ^^'^^ 

The last line is obviously invariant under a replacement of A and B hy CA and CB, 
respectively, where C G GL(2i?,C) is arbitrary. 

Now choose C{k) = {A + ikB)-^ G GL(2E, C) and find that (see also |KSn6bj ) 

C{k)A = -l{S{k) - 1) and C{k)B = {S{k) + l) . (3.8) 

Inserting this into (13.71) finally leads to 



S'{k) = ^{S{k) + t){Sik)-t) 

= -jj^isik) - siky) s{k) 



(3.9) 



which proves the statement for k eM.\ {0}. 

From Lemma [3. II we infer that S{k) is analytic in a neighbourhood of = and that 
the right-hand side of (13. 9p has a removable singularity at k = 0; hence (13.91) extends to all 
real k. Since, moreover, S{k) is unitary on R and analytic on C \ [±i(T(L)], and ^ as well 
as S{k)~^ are also analytic on this set, the full statement of the lemma follows by analytic 
continuation. □ 

From Lemma [3.11 we know that S{k) is meromorphic with finitely many poles on the 
imaginary axis. One can therefore perform power series expansions of the S-matrix in large 
parts of the complex fc-plane. Below we want to specify two such expansions, and to this 
end we exclude an annulus containing icr(L) \ {0}. This annulus is characterised by the 
two radii 

Amin/max := min/max{|A|; A G cr(L) \ {0}} . 

Here we assume that cr(L) ^ {0}; otherwise L = which is equivalent to the S-matrix 
being independent of k. 

We are now in a position to provide the announced expansions of the S-matrix. 
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Lemma 3.3. Let the same conditions as in Lemma \3. 1\ he given and assume that L ^ 0. 
Then the following power series expansions converge absolutely and uniformly in any closed 
subsets of the specified regions: 



1. For \k\ > A„ 



2. For\k\< Amin, 



oo 

S{A, B-k) = l-2P + 2Y,- {iir . (3.10) 



k 

n=l 



S{A,B;k) = -1 + 2P -2^2^^"" {^^y , (3.11 



n=l 



where P and L emerge from P and L, respectively, by replacing A, B with —B, A. 

Proof. For the first expansion we refer to the representation (13.51) of the S-matrix and 
employ the expansion 

X^-ik l + iXa . , oV^/'i^^A" 1 J 



Aq, -\- '\k 1 — tAo, 1 ^ ^ 



valid for \k\ > Aq,. Hence, for \k\ > Amax the S-matrix is 1 +2 Yl'^=i (iL)"- on rani?* and 
— 1 on keri?. Since L = on keri?, the relation f l3.10p follows. 

For the second expansion we remark that from (13. 3p one can readily deduce the relation 
S{A,B; k) = -S{-B,A; i) for k G M\ {0}, see |KS06aj . Thus, fl310D implies (KW\ . and 
the domain |A;| < Amin of convergence, in which S{k) is holomorphic, follows immediately. 

□ 

Lemma 13.31 also provides limiting expressions for the edge S-matrix as \k\ — > oo and 
I A; I 0, respectively, 

S^ = 1-2P and Sq = -1 + 2P , 

which we shall use subsequently. 

Later we shall integrate expressions containing the S-matrix along contours in the upper 
complex half plane and, therefore, we need to estimate the norm of the S-matrix along the 
contours. To this end we introduce 

_ /min{A G ct(L); A > 0} , if 3A„ > 
I OO , else 

and obtain the following. 

Lemma 3.4. Let A; G M and < k < A^^j^, then 

\\Sik + iK)l ||5(A:-iK)-i < max(l,^fii^| (3.13) 
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in the operator norm. Furthermore, if k > X 



max; 



then 



\\S{k + iK)\\, \\S{k-iK)-'\\ < 



K + X 



max 



(3.14) 



K — X 



max 



Proof. From (13.51) one can read off the eigenvalues of S{k + i/t) as ±1 and — (Aq, + k — 
^k)/{Xa — n + ik), a = 1, . . . ,d. In absolute values the latter quantities, as functions of 
/c G M, are maximised at /c = 0. Now suppose that A^;^ > and let < k < A^j„. Then 
the largest quantity among the | Aq, + k|/| Aq, — k| is the one on the right-hand side of (13.131) . 
In the case A^j^ = oo the upper bound is one. The proof for S{k — in)'^ is completely 
analogous. 

li K > Amax the same argument yields the bound (I3.14p . □ 

4 The spectrum of the Laplacian 

The scattering approach to the quantisation of a finite, metric graph utilises a secular 
equation based on the edge S-matrix of the graph. Here we closely follow the original 
approach as developed by Kottos and Smilansky |KS99b] for the case of (generalised) 
Kirchhoff boundary conditions, which was later generalised by Kostrykin and Schrader 
[KSOGbj . To keep the presentation sufficiently self-contained, we reproduce the relevant 
results below. We begin, however, with some general properties of Laplace spectra and 
finish this section with some remarks on the eigenvalue zero. 

4.1 Preliminaries on the spectrum 

Given a Laplacian on a compact metric graph, one would naturally expect that its spectrum 
is discrete and has a finite lowest eigenvalue. Kuchment indeed proved |Kuc04j that any 
such self adjoint (negative) Laplacian is bounded from below, and that its resolvent is trace 
class. Thus, the spectrum is discrete and bounded from below. Subsequently, Kostrykin 
and Schrader |KS06bj improved the lower bound. They showed that 



(4.1) 



where s > is the unique solution of 




(4.2) 



with /min denoting the shortest edge length and 




max{A G (r{L); A > 0} , if 3Xa > 
0, else 



(4.3) 
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Remark 4.1. As an aside we should like to mention that the lower bound (14.1 p is optimal 
in the following sense: Consider a trivial example of a metric graph given by an interval I 
of length /, and a Laplacian with domain specified by the choice 

A = AI2 and 5 = I2 , 

where A > 0, such that L = AI2 and o"(L) = {A}. (Equivalently, P = 0, Q = I2 and 
L = AI2.) Hence, (12.71) leads to the Robin boundary conditions 

H/(oJ + v-/'(oJ" ' 

and this implies the quantisation condition 

The solution k = in, with /t > A > 0, representing the lowest Laplace eigenvalue —k^ 
corresponds to a solution of the equation 

Ktanh — = A . 

This condition is equivalent to (14. 2p . demonstrating that the bound (14. ip is sharp for this 
'quantum graph'. 

Kostrykin and Schrader also showed |KS06b] that the number of negative Laplace 
eigenvalues is bounded by the number of positive eigenvalues of L (counted with their 
respective multiplicities). In the example above, the Robin Laplacian on an interval hence 
has at least one and at most two negative eigenvalues. 

For the trace formula one requires an a priori estimate on the number of eigenvalues. 
This is well known for the Dirichlet and the Neumann Laplacian, and in the case of Kirch- 
hoff boundary conditions can be found in |Sol02] . The same asymptotic law, however, 
holds also in the general case. 

Proposition 4.2. Given a self adjoint realisation of the Laplacian on a compact metric 
graph, the number of its eigenvalues G M counted with their multiplicities ) fulfils the 
following asymptotic law, 

N{K) ■=4^\j; <KH r^-K , K^oo, (4.5) 

where £ := /gi + ■ ■ ■ + is the total length of the graph. In particular, the Laplacian has 
infinitely many eigenvalues that only accumulate at infinity. 
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Proof. We prove the asymptotic law employing a variational characterisation of the eigen- 
values based on the quadratic form fl2.9p , as well as an analogue of the Dirichlet-Neumann 
bracketing (see |RS78j ). To this end we introduce two comparison operators. 

The first comparison operator is the direct sum of the Dirichlet operators on the edges, 
i.e., the Dirichlet-Laplacian of section [221 The domain of the associated quadratic form is 
characterised by the condition Fj,^ = 0, and therefore is contained in the domain of (12.91) . 
Moreover, on the Dirichlet-form domain both quadratic forms coincide. The comparison 
lemma devised in [RS78j hence implies that No{K) < N{K). Here N£,{K) is the counting 
function for the eigenvalues of the Dirichlet-Laplacian, which trivially fulfils the asymptotic 



As our second comparison operator we choose the direct sum of Robin Laplacians on 
the edges (see the example above). This operator is characterised by Pr = 0, Qr = I2E 
and L = XI2E with some A > 0. (A = in fact corresponds to Neumann Laplacians 
on the edges.) The boundary conditions are XFby + F^^ = 0, and thus they decouple the 
edges. The respective eigenvalue counting function Nfi{K) can be determined from (14. 4p 
and clearly obeys the asymptotic law (14.51) . The associated form domain, characterised by 
the condition PrF^^ = 0, contains the domain of (12. 9p . Choosing A = A^^^^, see (14. 3p . on 
the form domain of (12. 9p one finds Qr[F] < Qa[-^]- Therefore, the comparison lemma of 
jRS78] implies Nr{K) > N{K). 

Thus, NoiK) < N{K) < Nr{K) and the upper and the lower bounds both fulfil the 
same asymptotic law, which proves (14. 5p . The further statement follows immediately. □ 

This result is independent of any details of the quantum graph, apart from its volume. 
This is analogous to the corresponding results for Laplacians on manifolds or domains. In 
general, the asymptotic growth of the number of eigenvalues is proportional to the volume 
of the manifold/domain and to K^, where D is the dimension of the manifold/domain. 
This type of results is often referred to as 'Weyl's law' |Weyll| , and insofar Proposition 14.21 
is the quantum graph version of Weyl's law. 

4.2 The secular equation 

Apart from the edge S-matrix, the scattering approach requires the metric information of 
the graph, which enters through 



law (1431) . 




with t{l; k) : 



\ 




(4.6) 



where /c G C. Both matrices are then used to introduce 



U{k) := S{A,B-k)T{l-k) . 



(4.7) 



The topological and the metric data entering U{k) are hence clearly separated. 
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For real k the endomorphisms S{k), T{k) and U{k) of C^^ are obviously unitary. We 
therefore denote the eigenvalues of U{k) by e'^^^^^^ . . . , e^^^isCfc) Following Lemma [3.11 we 
conclude that U{k) can be extended into the complex /c-plane as a meromorphic function 
with poles at icr(L) \ {0}. The determinant function 



on which the scattering approach is based (see |KS99b] ). hence is also meromorphic. Its 
poles are in icr(L) \ {0}, but do not necessarily exhaust the entire set. 

Proposition 4.3 (Kostrykin, Schrader |KS06b] ). The determinant function (14. 8 p is mero- 
morphic on the complex plane with poles in the set icT(L) \ {0}. Furthermore, let kn G 
C \ [io"(L) U {0}] with ImA;„ > 0, then k"^ is an eigenvalue of —A, iff kn is a zero of 
the function (14.81) . i.e., F{kn) = 0. Moreover, the spectral multiplicity Qn of the Laplace 
eigenvalue k^ > coincides with the multiplicity of the eigenvalue one ofU{kn). 

Proposition 113] establishes a close connection between zeros of the determinant function 
(14. 8 p and Laplace eigenvalues. Notice that although Laplace eigenvalues occur as squares, 
k"^, the function (14.80 is not invariant under a change of sign in its argument. There exists, 
however, a functional equation under the substitution k i-^ —k. 

Lemma 4.4. For all C \ [±ia{L) \ {0}] the following identity holds: 



where M = E + d + dimker B and C = le^ + ■ ■ ■ + is the sum of all edge lengths. 
Proof. We decompose 



so that after taking determinants, 

F{k) = (-l)2^detT(A;) det S{k)F{-k) . 

Using the definition (14.60 of T{k) and the representation (13. 5p of S{k) then yields (14. 9p . □ 

We remark that when fc^ G M is a Laplace eigenvalue, then either G M when the 
eigenvalue is non-negative or, in the case of a negative eigenvalue, k = ii/t with k, > 0. 
The relevant zeros of the determinant (14.80 are then ±k G M or in, respectively. Unless 
K G <j{L) or /c = 0, all Laplace eigenvalues are covered by Proposition 14.31 However, in 
order to count Laplace eigenvalues in terms of zeros of F{k) with their correct multiplicities 
one has to establish a connection between the order of the zero and the multiplicity of one 
as an eigenvalue of U{k). In the trace formula we shall need this connection for the non- 
negative eigenvalues and hence now consider the eigenphases d{k) of U{k), defined through 



F{k) := det(l - U{k)) , 



(4.8) 




(4.9) 



1 - S{k)T{k) = -S{k) [1 - S{-k)T{-k)]T{k) , 



U{k)v{k) = e'^^^h{k) with \\v{k)\\ = 1 , 



(4.10) 
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for k eR. We then recall that U{k) = S{A,B; k)T{l- k) is analytic in C \ [\a{L) \ {0}]. 
According to analytic perturbation theory (see, e.g., |Kat95j ) its eigenvalues, for which we 
keep the notation e'^*-*^-', are continuous on this set and differentiable apart from possibly 
isolated points. Since, however, U{k) is real-analytic and normal for all A; G M, we can 
apply a sharpened version of analytic perturbation theory (see |LT85j ) to conclude that 
the eigenvalues are real analytic for all /c G M and that there exists a choice of eigenvectors 
with the same property. 

The following statement is a generalisation of a result found in |KS99bi IBWOSj that is 
valid for /c-independent S-matrices. 

Lemma 4.5. Let 6{k) be an eigenphase ofU{k), G M, with associated normalised eigen- 
vector v{k) = (t>i(A;), . . . , V2E{k)Y . Then 



2E 



dk 



1=1 



L^ + k^ 



v{k) 



Proof. Taking Lemma 13.21 into account we first observe that 



(4.11) 



U\k) = -^{S{k) - S*{k))S{k)T{k) + iSik)T{k)D{l) 



2k 
-2i- 



L 



where 



L^ + k^ 
Diil) 



U{k) + iU{k)D{l) 









D,{1) 



\ 



IeJ 



We also employed the relation 



2k 

that follows from (13. 4p . This then yields 



- Us{k) - S\k)) = -2\ ^ 



v{k),-[U{k)v{k)] 



-2ie'^^^^ (v{k) 



L^ + k^ 
+ {U*{k)v{k),v\k)) 



(4.12) 



v{k) ) + i{U*{k)v{k),D{l)v{k)) 



(4.13) 



(4.14) 



(4.15) 



-2ie'''('=) (v{k) 



L 



+ e 



L2 + fc2 

''^''^ {v{k),v'{k)) . 



v{k)) + ie''^^^ {v{k),D{l)v{k)) 



On the other hand, taking the derivative on the right-hand side of fl4.10p and then multi- 
plying with v{k) leads to 



ie'{k)e''^^^ + e''^^^vik),v'ik)) . 



(4.16) 
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Comparing fl4.15p and (14.161) proves the statement (14.111) . □ 

This lemma allows us to obtain an upper and a lower bound for the derivative of an 
eigenphase. Using /max/min to denote the largest and the smallest edge length, respectively, 
and introducing 



'^min ■ 



min{|A|; A G a{L) f] M"}, if BX^ < 
cxD, else 



in analogy to (I3.12p . we immediately get the following. 

Corollary 4.6. The derivative 9'{k) of an eigenphase 6{k) is bounded from above and 
below according to 

Un - TT- ^ ^'(^) ^ ^max + • (4.17) 

min mill 

In particular, if I min > '^/^min derivatives of all eigenphases are always positive. 
Proof. Obviously, 

2E 



1=1 



since the eigenvector is supposed to be normalised. Moreover, after a diagonalisation of L, 
when WLW* is diagonal with the eigenvalues A^ on the diagonal, one obtains 



r \ \ 



where w{k) = Wv{k). This first yields 



+ ' 



■^min \ + k^ I Aj^j^ 

and then finally (l4l7ll . □ 

As an important consequence of this corollary we are now able to extend the statement 
of Proposition 14.31 as required to count eigenvalues in terms of zeros of the determinant 
function. 

Proposition 4.7. Let the metric structure of the graph be such that l^vci > "^/^min- Then 
kf^ > is an eigenvalue of the Laplacian with multiplicity Qn, if ^kn G M \ {0} are zeros of 
the function F{k) of order Qn. 

Proof. From Proposition 14. 31 and Lemma S3] we know that the positive Laplace eigenvalues 
are in one-to-one correspondence with (pairs of) k values where U{k) has an eigenvalue 
one, and the respective multiplicities coincide. It hence remains to establish that the order 
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of the corresponding zeros of F{k) are exactly these multiphcities: From the definition 
(14. 8p of the function F{k) one obtains 



2E 

i=i 

so Corollary 14.61 imphes ^ (l — e^^^'-'^^) = —i9j(k)e^^^^''^ ^ 0, and the claim follows immedi- 
ately. □ 

We remark that since the scattering approach to the proof of the trace formula is 
based on counting zeros of the function F{k) on the real line with their multiplicities, the 
requirement /min > '^/^min essential. Otherwise one might count Laplace eigenvalues 
with incorrect multiplicities. Whenever L has no positive part, however, the condition is 
empty. This is, e.g., the case for non- Robin boundary conditions. 



4.3 The eigenvalue zero 

In general, zero is a Laplace eigenvalue as well as a zero of the determinant (14. 8p . and in so 
far Proposition 14.31 also applies to fco = 0. The spectral multiplicity go, however, typically 
is different from the degree of /cq = as a zero of F{k). For Kirchhoff boundary conditions 
it has been shown in jKN05] that the degree of the zero is E — V + 2, whereas the zero 
Laplace eigenvalue is non-degenerate, i.e., go = 1. Kurasov [KurOSj subsequently linked 
this difference in the multiplicities to the topology of the graph by noticing that a suitable 
trace formula contains the quantity 

and hence the Euler characteristic of the graph. This observation was generalised to yield 
an index theorem for any quantum graph with non-Robin boundary conditions by Fulling, 
Kuchment and Wilson |FKW07] . One can view the t-independent term ^ti S = ^{V — E) 
in the trace formula for the heat kernel (due to |Rot83] for Kirchhoff boundary conditions 
and |KPS07j for general non- Robin conditions) as a predecessor of this result. See also 
[BEOSj for a more detailed discussion. 

We here wish to give a further characterisation of the spectral multiplicity of the zero 
eigenvalue in the case of a general self adjoint realisation of the Laplacian. To this end we 



18 



first introduce, for /c G M \ {0}, the matrix 



C(/;fc) := 



/ U 



\ 



\ 



Ie 
2j:+lE 



H 
2i+lE 



2i+'i 



Ie 



(4.18) 



in which all matrix entries not indicated are zero. This now enables us to formulate the 
following. 

Proposition 4.8. For any given self adjoint realisation of the Laplacian specified through 
A,B, zero is a Laplace eigenvalue, iff one is an eigenvalue of S{A, B; k)C{i, k) for one, 
and hence any. A; G M \ {0}. Moreover, the multiplicity of this eigenvalue one coincides 
with the spectral multiplicity Qq of the zero Laplace eigenvalue. 

Proof. Eigenfunctions of the Laplacian corresponding to the eigenvalue zero must be of 
the form F = (/i, . . . , fs)^ with 



fj{x) = aj + PjX , X G [0, Ij] . 
Hence, the boundary values (12. 2p take the form 

Fbv = (ttl, • • • , CtE, ai + ...,«£; + Pih)'^ , 

Fi={P,,...,PE,-Pl,...,-PE)'^ . 

We now employ the boundary conditions (12.31) by using the expressions 
choice of A and B. The result can be rearranged to yield 

S{k)C4l;k) (^^^ =C4l;k) (^^ 

where, for any A; G M \ {0}, we have introduced 

with Di{l) as defined in (14.131) . We also use the abbreviations 



(4.19) 



for a possible 
(4.20) 



and f3 



\oieJ 



19 



The matrices C±{1, k) are invertible for all A; G M \ {0}, with 



C±{l;k)-' 



±1 
±2i-'i 



-1 



±2i-«i 



±2i~h/ 



We now substitute 



in (l4.2Up and obtain 



v{k) ■=C_{l-k) 



(4.21) 



S{k)C+{l]k)C^{l]k)-^v{k) =v{k) . 

It is straight forward to check that C+(Z; A;)C_(Z; A;)~^ = C{i,k), compare fl4.18p . The 
linearly independent eigenvectors of SC corresponding to the eigenvalue one then yield, 
via fl4.2ip and f l4.19p . coefficients aj and j3j for as many linearly independent Laplace 
eigenfunctions in L'^{T). □ 



We remark that the order of fco = as a zero of the function (14. 8p is the multiplicity 
of the eigenvalue one of 

U{0) - 







which is in no obvious way related to the multiplicity of the eigenvalue one of S{k)C{l; k) 
that appears in Proposition 14. 81 In the case of non- Robin boundary conditions, where the 
edge S-matrix is independent of k, however. Fulling, Kuchment and Wilson )FKW07j were 
able to relate the different multiplicities in the form of an index theorem. They showed, in 
particular, that then 

go--N=-tTS. 
As mentioned above, this term will reappear in the trace formula. 



5 The trace formula 

A trace formula expresses counting functions of Laplace eigenvalues in terms of sums 
over periodic orbits. Ideally, one would like to count Laplace eigenvalues fc^, with their 
multiplicities Qn, in intervals I in the form 

TTXi{-A) = J29n. (5.1) 
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The sharp cut-off provided by the characteristic function xi of the interval J, however, 
cannot be dealt with. One therefore replaces fl5.ll) with a smooth cut-off and, moreover, 
performs this count in terms of the associated wave numbers A;„, i.e., one seeks to find a 
representation for 

Y,9nKkn) (5.2) 

n 

in terms of sums over periodic orbits. One ambition then is to find a sufficiently large class 
of test functions h. It turns out that the following one parameter family of test functions 
is particularly suited for these purposes. 

Definition 5.1. For each r > the space consists of all functions /i : C ^ C satisfying 
the following conditions: 

• h is even, i.e., h{k) = h{—k). 

• For each h G H,. there exists 6 > such that h is analytic in the strip Mr+s '■= {k E 
C; I Imfcl < r + 5}. 

• For each h E Hj. there exists ?7 > such that h{k) = O (^■(J^^yT7;^ on Mr+s- 

We stress that the trace formula will only take non-negative Laplace eigenvalues into 
account. Hence, the test functions in (15. 2p will be evaluated at real arguments. We can 
therefore arrange the wave numbers /c„ > corresponding to Laplace eigenvalues k"^ in 
ascending order. Proposition 14.21 then readily implies the existence of a constant C > 
such that 

V gn\h{kn)\<CK sup \h{k)\ . (5.3) 
o<fc„<x '^elo.K] 

The third condition of Definition 15.11 now ensures that the sum in (15. 2p converges absolutely 
when h E Hr with some r > 0. Since this sum constitutes one side of the trace formula, 
choosing test functions h from any of the spaces Hr is an appropriate prerequisite for the 
trace formula. 



5.1 A precise eigenvalue count and the trace formula 

The trace formula is based on an exact count of non-negative Laplace eigenvalues with, 
however, a smooth cut-off as in fl5.2p . Essentially, this sum is expressed in terms of a sum 
over periodic orbits on the graph. In order to perform such an eigenvalue count one relies 
on the connection between Laplace eigenvalues and zeros of the determinant function (14.81) . 
More specifically, one chooses a test function h E Hr with some r > as well as some e > 0, 
which must be sufficiently small so that in the set C^.k = {k E C; | Im A;| < £, | Re k\ < K} 
the determinant function F{k) has only (finitely many) real zeros related to non-negative 
Laplace eigenvalues (compare Proposition 14. 3p . Then, by the argument principle, 

^{k)h{k)dk = Nh{0) + 2 9nKkn). (5.4) 

0<k„<K 
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Here N and g.^ are the multiplicities of /cq = and /c„ > 0, respectively, as a zeros of F{k). 
Following Proposition 14.71 gn also is the multiplicity of the Laplace eigenvalue fc^. The 
factor of two occurs since kn, —kn G Ce^k and we have exploited the fact that h{k) is even. 
Based on this relation the trace formula emerges when one expresses F'{k)/F{k) in terms 
of a suitable series (eventually leading to a sum over periodic orbits), and performs the 
limit K oo. The result of this procedure is summarised in the following statement. 



Proposition 5.2. Let /mm > '^/■^mm '^'^'^ choose h G Hj. with any r > 0. Then 
Nh{0) + 2Y,9nhik^) = J2^J tT[Aik)U\k)] h{k) dk , 



(5.5) 



n=l 

where 



O T 



Proof. Our strategy is to show that both sides of fl5.5p equal 



lim 

£^0+ 27ri 



+ 00 



F' F' 

— {k — ie) h{k — ie) {k + ie) h{k + ie) 

F F 



dk . (5.7) 



Beginning with the left-hand side, we have to show that in 05.41) the limit K ^ oo can 
be taken, followed by e — > 0. To this end we notice that the right-hand side of 05.41) is 
explicitly independent of e, in the range described above that equation. Furthermore, from 
05. 3p we already know that the sum over kn converges absolutely in the limit K ^ oo. 
In order to perform these limits on the left-hand side of 05. 4p . and thus producing 05.71) . 
we have to estimate the contribution to the integral coming from the vertical parts of the 
contour. 

Lemma F3.3I implies that for G C with \k\ > Amax the approximation F{k) = Foo{k) + 
0{\k\~^) holds, where 

F^{k) := det{l - SooTik)) . (5.8) 

Since F^o depends on k only through the matrix entries e''^'^ of T{k), see 04. 6p . it can be 
represented as 

b 

Foo(A;) = l + ^rf„e'^"\ (5.9) 

n=l 

with some 6 < oo. Moreover, /3„ > is a finite sum of edge lengths, and (i„ G C is an 
appropriate coefficient. The expression 05.90 is first defined for A; G M, but can be readily 
extended to complex k. Since S'{k) = 0{\k\~^), see Lemma \^?2\ and eq. 04.140 . we also 
find F'{k) = F^{k) + 0{\k\~^). For sufficiently large \k\ one can therefore approximate 
F'{k)/F{k) by F^{k)/ Foo{k). In order to estimate the latter we employ 05.90 to obtain 

|i^^(fc)| < &rfmax/5maxe-^-'°^'"'^ , (5.10) 

with (imax := niax{|(i„|} and /5max/min := max / min{/?„}. Note that this bound is indepen- 
dent of Rek. Furthermore, we pick fc(o) G M such that F(A;(o)) 7^ as well as -Foo(^(o)) 7^ 0, 
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and take advantage of the fact that Foo{k), A; G M, is an almost periodic function (see, e.g., 
[Boh47] ). One can hence construct a (strictly increasing) sequence {k^y, j G Nq} with 

|Foo(%)) - F^(fc(o))| < , J > . (5.11) 

Hence, in particular, Foo{k(^j)) 7^ 0. The estimate fIS.lOp . moreover, implies that 

|Foo(%)) - Foo(A:o-) + < (5.12) 

when \k,\ is sufficiently small. Therefore, when e is small enough, the function |Foo(fc(j)+ifi;)|, 
—e < K < e, is uniformly bounded from below away from zero. Thus there exists Cs > 
such that 

< Ce snp{\h{k(^j) + iK)\; kE (-£, +e)} . (5.13) 

The third property of test functions required in Definition l5.1l then ensures that this integral 
vanishes as fc^) —>■ 00. A completely analogous reasoning applies when replacing k(^j^ 
with —kf^j). Hence, as i^' — > 00 the contributions coming from [—K + ie, —K — ie\ and 
[K — ie,K + ie\ to the integration along dC^^K in (15.41) vanish, and the limits e and 
K 00 can be interchanged. 

It remains to prove the equality of (15. 7p with the right-hand side of (15.51) . In order to 
achieve this we first recall from Proposition 14.31 that the function F{k) is holomorphic and 
zero-free in the strip {/c G C; < ImA; < ei}, where ei := minjA^j^, Ki} when —nl is the 
largest negative eigenvalue of —A; if —A > we set ei := Hence, 

^[k + ie) = ^logdet(l - U{k + ie)) = -tr(^(l - f/(A; + ie))"V'(A; + ie)) , 
when < 6 < El- Furthermore, following Lemma [3. 131 we obtain the bound 

\\U{k + ie)\\<\\S{k + ie)\\\\T{k + ie)\\<max\l,^^^^\ e^"'-'" (5.14) 

I ^min ~ ^ J 

in the operator norm. Thus, when 

e < A+„tanh(^^y^) , (5.15) 

the operator norm of U{k + ie) is less than one. Iff /min > 2/A^j„ this condition is solvable, 
providing some £2 with < £2 < Aj^j^ such that (I5.15P is true for all e G {0,62). For such 
e the expansion 

00 

(1 - U{k + ie)y^ = J2u{k + ieY 

1=0 
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holds. Moreover, a straight forward calculation based on the relation (14. 12p produces 
tT[U{kyU'{k)] = tT[A{k)U{ky^^], see f l5.6p . Therefore, after cyclic permutations under the 
trace, 

TP' °° 

— {k + ie) = -^tr(A(A; + i£)[/'(A; + i5)) (5.16) 
1=1 

is true for all G M and e G (0,63), where £3 := min{ei,e2}- 
Likewise, based on the relation 

(1 - U{k - ie)y^ = -(1 - U'\k - ie))~^U'\k - \e) 



that holds for < £ < Ai;^, we obtain 

00 

\k-\e) = ^tx{k{k-\e)U~\k-\e)) (5.17) 



min' 

F' 



1=0 



for A; G M and sufficiently small e > 0. 

We now want to use the representations 05.161) and (15.171) in (15. 7p . interchange integra- 
tion and summations, and finally perform the limit e ^ 0. In order to achieve this we first 
notice that (for fixed e) the estimate 

\tT{A{k±ie)U^\k±ie)) \ < 2E\\A{k ± ie)U'^\k ± ie)\\ 

< 2E\\A{k ± ie)\\\\U^\k ± ie)\\' ^^'^^^ 

justifies to interchange integration and summation according to the dominated convergence 
theorem. Next, the contour of the integral 



Im fc=it£ 



tT[A{k)U^^{k)] h{k) dk (5.19) 



can be deformed into Im A; = ±5 with a sufficiently small 5 > 0; in particular, 6 < r. We 
then fix 6 and find 

f tT[A{k)U^\k)] h{k) dk 

J Im k=±e 

/+00 
|tr(A(A; ± iS)U^^{k ± 16)) \ \h{k ± i6)\dk (5.2O) 
-00 

:^mm^g-5U. / \\A{k±i6)\\\h{k±i6)\dk , 
-^min ~ ^ J J -00 

when using (15.141) and (I5.18P with 6 instead of e. Since the conditions in Definition 15.11 
apply, the integral on the right-hand side is finite and the positive constant raised to 
the power / is smaller than one. The sum on / hence possesses an absolutely convergent 
majorant, uniform in e in the range indicated. Thus the summation and the limit e — can 
be interchanged. Furthermore, another application of the dominated convergence theorem 
allows to perform the limit £ of (I5.19P inside the integral, finally proving (15.51) . □ 
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In order to arrive at the trace formula itself, the sum on the right-hand side of fl5.5p 
has to be reformulated as a sum over periodic orbits. The summation index / then denotes 
the topological length of the orbits and the trace of K{k)U\k) is identified as a sum over 
the set Vi of periodic orbits with topological length I. The term with / = plays a special 
role and will be treated separately, whereas contributions with negative / are related to 
those with —I in a simple way. 

Before we state the trace formula, however, we introduce some notation: The 'volume' 
of the graph is the sum C = 1^^ + ■ ■ ■ + 1^,^ of all edge lengths. Furthermore, ii h & Hr is a. 
test function, its Fourier transform is 



1 r+oo 

h{x) = — / h{k) e''^^ dk . 
2vr J-oo 



We recall that the second and the third property required for h in Definition [STT] guarantee 
that h{x) = 0(e~™) as a; — > cxd. Moreover, the Fourier transform of a product A{k)h{k) is 
the convolution of the respective Fourier transforms, i.e., it reads 

/+00 
A{x -y)h{y) dy = A*h{x) . (5.21) 

Below this convolution will often have to be understood in a distributional sense, as the 
functions A{k), though being regular, not always decay sufficiently fast as — > oo. 
We are now in a position to state the first variant of the trace formula. 

Theorem 5.3. Let T he a compact, metric graph with a self adjoint realisation — A(y4, B\ I) 
of the Laplacian, such that the condition Ijam > is fulfilled. Furthermore, let h G Hj. 

be a test function with an arbitrary r > 0. Then the following identity holds: 



^ 1 r+oo 

J2 9n hikn) = C hiO) + {go - |iV) hiO) - — / hik) 

oo 



ImtTS(k) 
dk 



(5.22) 



1=1 P<^Vt 



Here Ap is the Fourier transform of the amplitude function Ap{k) associated with every 
periodic orbit p. This function is meromorphic with possible poles at the poles of the S- 
matrix, and has a Taylor expansion 

oo 

Mk) = ^afk-^ (5.23) 
j=o 

that converges for \k\> Amax- In general, only the sum over the topological lengths I on the 
right-hand side of (I5.22p converges absolutely, but not the entire double sum over periodic 
orbits. 
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Proof. We have to evaluate tr[A(A;)f/'(A;)], and first notice that according to fl4.6p and fl4.7p 
the matrix entries of U are 

where 

^ ^ fj + E, iil< j<E 

' -E, ii E + l< ] <2E ' 

Hence, the indices j and Uj correspond to the two edge ends of the j-th edge. Local 
boundary conditions then imply that Sj^^^^^ ^ requires the edges with ends ji and to 
be adjacent. Therefore, when / > the non vanishing terms in the multiple sum 

2E 

tr(Af/') = KnSn^,, ■ ■ ■ ^..^.o e'^^^'-^--^'-) (5.24) 

jO,..;jl=l 

correspond to the closed paths of topological length / on the graph. 

We then make use of the decomposition (15. 6p of A and begin with the contribution 
of tr[D(/)[/'(fc)], which can be evaluated as in the case of Kirchhoff boundary conditions 
|KS99b] : Due to the specific diagonal form of D the terms in (15.240 corresponding to closed 
paths related by cyclic permutations of their edges can be grouped together. This finally 
yields a sum over the periodic orbits of topological length /, 

tr[D(Of/'(fc)] =2 5^/#yli,p(A;)e''=''' . 

According to (15.241) . the functions Ai^p{k) result from multiplying the local S-matrices of 
the vertices visited along the periodic orbit p. Moreover, is the primitive length of p, 
i.e., the length of the primitive periodic orbit associated with p. Due to the relation (I5.2ip 
we therefore obtain 

1 r+oo 

— / tr[D(Z)[/'(A;)]/i(A;) dA; = 2 5^/#(/i*ii,p)(g . 

The case of negative I follows by noticing that 

U-\k) = T-\k)S~\k) = T{-k)S{-k) , 



and thus 



tr[D{l)U~\k)]=2j2ltAiAk)e-'''^ , 



leading to 

— / tr \D{l)(U\k) + U-^{k))] h{k) dk 



2^/#[(/i*4p)(g + (/i*Ai,p)(/p 

pep. 



26 



In order to calculate the contribution from tr'^ ^ 
L 



-2itr 



U 



ii[S'TU 



2E 
Jiv,i;=l 



jj^U'' (k)] we notice that 

mn+-+hi) . (5.25) 



Again, the multiple sum can be viewed as a sum over the closed paths of topological length 
I, and the contributions of representatives of periodic orbits can be grouped together. 
Eventually, this leads to 



- 2tr 



L2 + fc2 



U\k) 



(5.26) 



where the functions A2^p{k) emerge from multiplying local S-matrix elements and their 
derivatives along the closed paths as specified in fl5.25p . Negative I are dealt with as above, 
so that the contribution from I E'L \ {0} to the sum on the right-hand side of (15.51) yields 
the sum on the right-hand side of (I5.22p . with Ap{k) = lfAi^p{k) + A2,p{k). 
The contribution coming from / = finally is 



1 
27ri 



trA{k)h{k) dk 



1 

~2^ 
2Ch{0) 



+00 



tr 
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+ 00 



h{k) dk 



Adding the contribution of the Laplace-eigenvalue zero to (15. 5p . after a rearrangement of 
the terms the result (15.220 follows. □ 

As mentioned in Theorem 15.31 the double sum over periodic orbits in (I5.22p does not 
converge absolutely when h E Hr with arbitrarily small r > 0. In the following we are going 
to show that an absolutely convergent periodic orbit sum can be achieved under sharpened 
conditions on the test function. In order to formulate these conditions we introduce the 
function 

-log(2E) + -artanh - 



1{K) 



K 



VA 



< fi: < A+ ^ , 



(5.27) 



which attains its unique minimum at some a G (0, A^^j^). Moreover, the minimum can be 
bounded from below as l{a) > (2 + log(2i?))/Aj^i^ so that, in particular, l{a) > 2/A^jj^. 

Theorem 5.4. Let T be a compact, metric graph with a self adjoint realisation —A{A, B; I) 
of the Laplacian, such that the condition /min > ^(o") is fulfilled. Furthermore, let h G Hr 
he a test function with an arbitrary r > a. Then the following identity holds: 



CO „ 

J2 9n h{k^) = C h{0) + {go - iAT) h{0) - — / 

n=0 ^ 



— oo 



k 



+ ^y{h* Ap){ip) + {h*Ap){ip) 



(5.2J 
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The quantities appearing are the same as in Theorem \5.3[ Here, however, the sum over 
periodic orbits converges absolutely. 



Proof. The task is to analyse the convergence of the sum over periodic orbits in f l5.28p 
more closely. In Theorem 15.31 we only considered the convergence of the sum over /. 
This was based on the estimate f l5.20p . which shall now be refined. To this end we refer 
to (15.241) and fl5.25p . and notice that ti[K{k)U\k)\ is a sum consisting of (2i?)' terms 
Siiiiik) ■ ■ ■ Si^ii^^j^{k)h{k), each of which contains a product of / factors that are matrix 
elements of either S{k), iD{l), or ijj^- These factors are all holomorphic in a strip 
< Irak < A^jjj and, following Lemma [3. 4[ in any strip <lmk < k, with k, < A^j^, they 
are bounded from above in absolute value by a constant times (A^J^j^ + k)/ {X^^^ — n) . Thus, 
in particular, for all k < min{A^i^,r} any product Si^^i^i^k) . . . Si^i^_^-^{k) h{k) is holomorphic 
for < Im k < K. Since here we are evaluating Fourier transforms at /j^ + ■ ■ ■ + /j^ > 0, we 
can apply a suitable version of the Paley- Wiener theorem that is concerned with Fourier 
transforms of holomorphic functions in vertical strips of the upper half plane. This is an 
obvious, slight variation of Theorem IX. 14 in |RS75j . and implies that there exists > 
such that 



+00 2_E 

ii[K{k)U\k)\ h{k) dk < J2 



+00 

s,,,i,(fc)...s,,,,(fc)Mfc)e^'=('n+-+'0 dk 

■oo 

I 



(5.29) 



where we made use of the function defined in (15.271) . We hence conclude that under 
the conditions stated the sum over periodic orbits converges absolutely. In particular, 
^min > ^(^) must be fulfilled. The latter condition is optimised for k = cr, at which the 
function 1{k) attains its minimum. For the above to hold, the condition h & Hr with r > a 
must be satisfied. □ 

We remark that whenever an amplitude function Ap{k) associated with a periodic orbit 
is independent of k, as it is always true for non- Robin boundary conditions, the convolution 
h * Ap{lp) degenerates into a product Aph(lp). In any case, the condition r > a imposed 
on the test function h ensures that h * Ap{x) = 0{e~'^^^^), which enables the absolute 
convergence of the periodic orbit sum. 



5.2 The trace of the heat kernel 

The first trace formula of a quantum graph, due to Roth |Rot83j . expresses the trace of the 
heat kernel for the Laplacian with Kirchhoff boundary conditions in terms of a sum over 
periodic orbits. This has recently been extended to all non-Robin boundary conditions by 
Kostrykin, Potthoff and Schrader |KPS07j . An application of Theorem 15.41 now allows us 
to produce an appropriate trace formula for all self adjoint realisations of the Laplacian. In 
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such a case, however, negative Laplace eigenvalues —n^, with multiplicities , may occur 
so that the trace of the heat kernel, 

Tr e^* := 9n e'^"* + E ' ^ > ^ ' 

-kI<0 kl>0 

does not follow immediately from Theorems 15.31 and 15.41 Instead, we consider 

Tr+e^* := E ^^e-'^"* , t>0. 

In the trace formula we hence have to choose the (entire holomorphic) function h{k) = e~'^^*, 
t > 0, which is in Hr for all r > 0. Its Fourier transform is h(x) = ■ We also 

introduce the functions 

"+00 



1 

a,{t) := -= / A,[l, - y) e"^'/^* dy 
V47rt J-00 



associated with the periodic orbits on the graph and obtain the following statement. 

Theorem 5.5. Let T be a compact, metric graph with a self adjoint realisation —A{A, B; I) 
of the Laplacian, such that the condition /min > l{a) is fulfilled. Then the following identity 
holds: 

Tr+e^* = + {g, - ^N) - i^i^^e^"* erfc(|A,| v^) + 2 Re a,(t) . (5.30) 



* a=l ' pGV 

Here d is the number ofnon vanishing eigenvalues Aq, of L (counted with their multiplicities) 
and 



erfc(x) = —= / e ^ dy , x > 



is the error function complement. Moreover, as t ^ 0^ the trace of the heat kernel has a 
complete asymptotic expansion in powers of v^, whose leading terms read 

Tre^* = -^ + 7 + 0(Vt) , t 0+ . (5.31) 

Here 

N ^ _ 1 ^ Iv-^ d+-d^ 
7 = 9o - Y + 1^9n - 2 Z^^OJ + 2Z^^Pd 2 ' 

n j 3 

where and 'jpj are the orders of the finitely many, non zero, purely imaginary ze- 
ros and poles, respectively, of the determinant function F and d± denotes the number of 
positive/negative eigenvalues of L. 
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Proof. We first observe that the expression fl4.14l) implies 



2kL _ A A. 



Imtr^(A;)=tr-^^ = 2A:}_^- 



a=l 



A;2 



then we employ the representation 



,,2 



2 X _ 2 f e~y 
erfc(|a;|) = e"^ / — — — 5- , xG 



Jo + 

Hence the relation (15.301) is an immediate consequence of (15.281) . when the test function h 
is chosen as indicated above. 

In order to determine the small-t asymptotics we go back to the relation (15.41) . in which 
we use the test function h{k) = e~'^^*, t > 0. We then deform the contour into dCf^^K, with 
[3 > max{Amax, s}, where s from (14.21) is such that — yields a lower bound on the Laplace 
spectrum. Thus the contour now encloses all non real zeros and poles of the determinant 
function F and, therefore, in this process we pick up contributions from all poles of F' / F 
on the imaginary axis. 

Having to perform the limit K 00 with /3 kept fixed, we need to estimate the 
contribution coming from the vertical parts of the contour, i.e., for | Re/c| = K and e < 
|ImA;| < (3. Firstly, F{k) is a polynomial in the matrix entries of S{k) and T{k). The 
latter are e''^'*, whereas the fc-dependence of the former is given by (A^ — iA;)/(AQ, + ik), 
see (13.51) . In the two strips e < |Imfc| < /?, with neighbourhoods of the poles at iA^ 
removed, all matrix entries are bounded, and hence F'{k) is of polynomial growth in k. 
Secondly, for sufficiently large \k\ we again approximate F{k) by Foo{k), see (15.81) . and 
write Foo{k) = Iljfill ~ ^i(^))- Here ui{k), . . . ,U2E{k) are the eigenvalues of SooT{k), 
which satisfy 

\ujik)\ < \\S^\\ ||T(A;)||<e-^'-'"<l, 
for all k in the strips e < \ lmk\ < p. Thus, in these strips, 

\F^{k)\> (l-e-^'-'")'''>0. 

Thirdly, when \K\ > (3 the factor e~*'^^ is of the order e~*^^ in the strips. Hence, the 
integrand of (15. 4p on the vertical parts of the contour with | ImA;| > e is bounded by a 
polynomial times e~*^ . As K ^ 00, these parts of the contour therefore do not contribute 
to the integral. Hence, from (15.41) we obtain 



N + 2Y,9n e-^"* = ^v. e'^^'^* - ^ 7o. e^o 

n=l j j 



1 

27ri 



+CXD 



F F 



&k . 



Here the non real zeros of F are denoted as i/toj and its poles as i/tpj; the respective orders 
are 70 j and 7^^-. 
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To proceed further we follow the argument leading from fl5.14p to (15.161) and fl5.17p . as 
well as the subsequent discussion of interchanging integration and summation. As we are 
dealing with large /3 instead of small e, due to fl3.14p the estimate f l5.14p is replaced by 



||f/(A; + i/3)|| <max<'l,|^^^l^ e"^' 



mill 



Hence, for A; G M and (3 sufficiently large analogous relations to f l5.16p and fl5.17p are 
obtained, eventually leading to 



n=l j j 



+ ^ — / tr [A{k + i/3)U\k + 1/3)] e-(^+'^)'* dA; 
1=1 

+ y — / tr \A{k - ip)U-\k - iP)] e-^^-'P)"' dk . 

ZTTl / 

/ — n ^ — CXJ 



In the integrals with / 7^ we now replace /? by /?/v^, < t < 1, and change variables 
from k to q = k\/i, yielding 

1 r+00 

I^{t,f3) := ^ j ^ tr[A(^(g±i/5))[/±'(^(g±i/3))] e'^^^^^)^ dq . 
These integrals can be bounded in analogy to f l5.29p . 



Summing over / 7^ then finally shows that these contributions can be estimated as being 

The term with / = can be calculated explicitly and yields the same contribution to 
the heat trace as the first term and the sum over a{L) \ {0} on the right-hand side of 

( Km . 

We add the contributions of negative Laplace eigenvalues and use that erfc(a;) has a 
complete asymptotic expansion in x, with erfc(a;) = 1 + 0{x), as a; — >• 0. The expansion 
f l5.3ip then follows immediately. □ 

At this point we recall that in the case of non- Robin boundary conditions 7 = go — \N = 
^trS*, which has also been given an interpretation as (one half of) a suitable Fredholm 
index. In this case, therefore, the constant term in the small-t asymptotics of the heat 
kernel has a topological meaning. 

Finally, we should like to mention that a suitable Tauberian theorem (see, e.g., |Kar31j ) 
allows us to recover Weyl's law 

N{K) r^-K, K ^00 , 

TX 

see also Proposition 14. 2[ from the leading term in the expansion (15.310 . 
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6 Conclusions 



Our principal goal was to investigate spectra of general self adjoint realisations of Laplace 
operators on compact metric graphs, culminating in proofs of some trace formulae. In this 
context we achieved to allow for a large class of test functions, leading either to absolutely or 
conditionally convergent sums over periodic orbits, respectively, representing appropriate 
spectral functions of the Laplacian. 

Previous work on quantum graph trace formulae |Rot83l IKS99bl IKN051 IKPSOTj was 
restricted to Laplacians with non-Robin boundary conditions. As compared to these cases 
there are some modifications we had to take care of. Firstly, non-Robin boundary condi- 
tions correspond to /c-independent S-matrices and hence do not involve any derivatives of 
S{k). This is in line with the fact that L = 0, so that in the trace formula the contribu- 
tion (15.261) to the periodic orbit sum is absent. Moreover, since therefore \^^^ = oo, the 
restriction imposed on l^i^ in the general case is void so that any set of lengths can be 
attributed to the edges. What still remains to ensure an absolutely convergent periodic 
orbit sum in the case of non-Robin boundary conditions is the single requirement h & 
with r > (log 2£') //mill on the test functions. This implies h{x) = 0(e~'''^'), which in turn 
compensates for the growth in the number of periodic orbits entering the sum 

p£V,ip<e 

when £ — > oo. 

Another property of Laplacians with non-Robin boundary conditions is that they are 
non-negative. This fact is linked to the non-positivity of L which, as L = 0, is trivial. The 
determinant function (14.81) hence is entire holomorphic with only real zeros in the complex 
non-negative half plane, see Proposition 14.31 

As we have shown in Theorem 15.31 the condition on the test functions can be relaxed 
to h E Hr with any r > 0, so that h{x) = 0{e~^^^^) with some (arbitrarily small) 5 > 0, 
when one is willing to accept a conditionally convergent sum. This has to be understood 
in the sense given in (I5.22p . i.e., where the terms are arranged as a double sum over the 
topological lengths of the orbits and over the periodic orbits of fixed topological length. 
For non-Robin boundary conditions a refined analysis of convergence had produced even 
more relaxed conditions to be demanded from the test functions, see |Win08] . 

An important application of the trace formula for quantum graphs with non-Robin 
boundary conditions was to prove an inverse theorem, very much in the sense of Kac's 
famous question 'Can one hear the shape of a drum?' [Kac66j . Gutkin and Smilansky 
[GSOl] showed that under certain conditions, which include the requirement that the edge 
lengths be rationally independent, the Laplace spectrum determines the connectedness and 
the metric structure of a compact metric graph uniquely. In this sense isospectral quantum 
graphs are isomorphic. Gutkin and Smilansky made essential use of the trace of the wave 
group, 

oo 

Tr e""^ + c.c. = 2 ^ ^„ cos(A;„t) , t , 

n=0 
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which can, in the case of non-Robin boundary conditions, be expressed as a sum of 6- 
singularities at the lengths Ip of periodic orbits. Searching for these singularities then 
allows one to first identify the geometric length spectrum. In a second step a certain 
algorithm can be used to determine the connectedness and all individual edge lengths. 
After a slight modification this proof can now be taken over to the case of Robin boundary 
conditions almost verbatim. For this purpose one reads off from the trace formula fl5.28p 
the distributional identity 

CO _ d , 

J2 9n CO^{Kt) = C 6{t) +g,-lN--J2TYl 

„ ^ -I \^a\ 

n=u a=l 

+ Re[ip(/p -t) + Ap{lp + t)] . 

Since, in general, the periodic orbit amplitudes are functions of k, there are no longer 
pure (5-singularities present at the lengths of periodic orbits. For large k the amplitudes, 
however, possess the expansions fl5.23l) with leading terms These do not vanish since 
they stem from the corresponding leading term S^o = 1 — 2P of the S-matrix. Hence, the 
Fourier transforms Ap of the amplitudes have leading singularities of 5-type at the lengths 
Ip of periodic orbits. Therefore, applying the algorithm of Gutkin and Smilanksy to this 
wave trace enables one to identify the graph connectivity and the edge lengths in the same 
way as previously. 

In summary, almost all spectral properties established so far for Laplacians on compact 
metric graphs with non-Robin boundary conditions carry over to arbitrary self adjoint 
realisations of the Laplacian. Therefore, there are many more quantum graphs model 
available that are suitable for further investigations in, e.g., the field of quantum chaos. 
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